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Learning Numerical Analysis through 
DERIVE T. Etchells, J. Berry 

This book covers the major numerical methods, and 
their analysis, for first courses at college and under- 
graduate level. The relative merits of each method 
are covered both analytically, providing a thorough 
grounding in the algebraic approach, and practically, 
through the tried and tested computer lab-based ac- 
tivities. 

DERIVE provides a platform on 
which to quickly and accurately 
perform many complicated numeri- 
cal calculations. Also, DERIVE’s 
ability to algebraically manipulate 
expressions and perform calculus 
operations, enhances the investiga- 
tion of the convergence of numeri- 
cal methods. Each chapter includes 
the development and algebraic 
analysis of the methods, lab-based activities, ideas for 
coursework, case studies, exercises and solutions. 
Free supporting utility files are downloadable via 
Chartwell-Bratf s web server. 

Chapter 1 introduces the basic tool of numerical 
methods, which is recurrence relations, their solution 
and ill-conditioning problems. In chapter 2 we use 
recurrence relations methods that are used in solving 
equations. Chapter 3 deals with the approximation of 
functions by polynomials, and in particular the Taylor 
Polynomial, which is then used extensively in chapter 
4 to analyse the errors associated with numerical 
methods. Chapters 5 and 6 deal with numerical ap- 
proaches to the calculus of differentiation and inte- 
gration. In chapter 7 we introduce and analyse nu- 
merical methods of solving differential equations. 
ISBN 0-86238-468-0, 239 pages, Chartwell-Bratt, 
February, 1997. 


Mathematical Activities with Computer 
Algebra, a photocopiable resource book from 
Etchells T, Hunter M, Monaghan J, Pozzi S and 
Rothery A 

This photocopiable resource book is the first of a new 
generation of support materials for the educational 
use of computer algebra. Designed to be used with 
any computer algebra system, the authors go beyond 
mere button pressing and show how 
to harness the power of computer 
algebra systems for educational 
purposes. 

Concepts are illustrated, techniques 
and methods presented, and model- 
ling and applications are explained. 
Appendices give overviews of 
DERIVE, Maple, Mathematica, 
Theorist (MathPlus) and the new TI- 
92 calculator. 

Activity Worksheets, Help Sheets 
and Teaching Notes cover a wide range of mathe- 
matical topics at school and college level. 

Topics covered include; functions and graphs, differ- 
entiation, integration, sequences and series, vectors 
and matrices, mechanics, trigonometry, numerical 
methods. Activities include: Multiplying factors; 
Equation of a tangent; Taxing functions; The tile 
factory; Function and derivative - visualisation; The 
approximate derivative function; Sketching graphs; 
Pollution and population; Max cone; Optimising 
transport costs; Area under a curve; Enclosed areas; 
A function whose derivative is itself; Wine glass 
design; The limit of a sequence; Visualising Taylor 
approximations; Visualising matrix transformations; 
Blood groups; Circular motion; Swing safety; No 
turning back; Modelling the sine function; Solving 
equations with tangents. 

20 Pounds, 96 pp, 1996, ISBN 0-86238-405-2 



Both books are available in 201 1 (e.g. AMAZON), Josef 

In the last DNL I gave a wrong ISBN for 'An Introduction to the Mathematics of Biology'. Here is the 
correct one: ISBN 3-7643-3809-1 . 


Interesting WEB sites 

http://www.swp.co.at/swhe/swhe.html (Soft Warehouse Europe) 
http://derive.com (Soft Warehouse Hawaii) 

http://ti.com.calc (for TI-92 programs); http://www.studli.se/chartwell.html (Publisher) 

http://www.kolleg.nuernberg.de/ti92.htm (W.Propper's Tl-programs: ableit(), galton(), binvert()) 
http://www.tech.plym.ac.uk/maths/CTMHOME/CTM.HTML (John Berry in Plymouth) 
http://www.cms.livjm.ac.uk/www/homepage/cmstetch/index.htm (Terence Etchells, DERIVE and Tl) 
All these sites are now (201 1 ) only of historical interest, Josef 

Do you know other sites? Share your knowledge with us! We will visit some DUG members' home 
pages in the future. Jan Vermeylen sent a valuable collection of addresses. Wait for the next issue!! 
Thanks to Jan from Kapellen in Belgium. Here is one of his goodies: 

http:Warchives.math.utk.edu (Mathematical Archives with among others a lot of historical information) 
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Liebe DERIVE- und Tl-Freunde, 

Sie werden sicher uber das Bild erstaunt sein. 
Ich bin in keinen Jungbrunnen gefallen. Das ist eines 
der jungsten Mitglieder der DERIVE-Gemeinde: Kim- 
berly, unsere erste Enkelin im zarten Alter von einer 
Stunde. Damit sind nun we/e Fragen beantwortet 
(Hallo, Barbel & Co). 

Hier ist nun unsere 25. Ausgabe und ich habe 
aus AnlafS dieses kleinen Jubi- 
laums SWHH und SWHE um 
einen Beitrag fur den DNL#25 
gebeten. Es sind zwei Gruflworte 
geworden (Seite 37), fur die ich 
mich herzlich bedanke. Ich dru- 
cke sie gerne ab, mochte aber 
die vi e/e Anerkennung, die darin 
ausgedruckt wird, gerne mit Ih- 
nen alien teilen. 

Manchmal kommen aller- 
dings auch leise Klagen: .... die 
Themen - und die zugehorigen 
Losungen - sind hin und wieder 
von einem anderen Stern, sehr 
weit weg von meinen taglichen 
" Problemchen " mit der Mathema- 
tik (D.BIum) 

Dazu mochte ich bemerken, 
dad ja Gott sei Dank nicht nur Lehrer zu unseren 
Mitgliedern zahlen und dad es auch den Lehrern ges- 
tattet sein sollte so dann und wann nach den Sternen 
zu greifen. Ich mochte aber alien, die sich auch betrof- 
fen fuhlen, ermuntern, uber Ihre "Problemchen" zu 
berichten. Da kann keine Zuschrift und kein Beitrag zu 
wenig aufregend Oder zu einfach sein. Sie alle gestal- 
ten den DNL. Und glauben Sie mir bitte, auch ich 
habe meine Mathe-"Problemchen" mit den 15- 
19jahrigen, und das oft nicht zu knapp. 

Ich habe mich aber bemuht, gerade diesen DNL 
besonders fur den Lehreralltag nutzbar zu machen. Im 
nachsten DNL werden wir uns auch der Geometrie 
widmen: Abbildungen von Objekten im R 3 , mit ver- 
deckten Kanten, Fraktale usw. Auch eine ganze 
Schachtel voll mit Tl- Fragen und Antworten steht fur 
Sie bereit. 

In vielen Fallen konnen DERIVE-Beitrage An re- 
gunge n fur den Einsatz mit dem Tl geben und umge- 
kehrt. Wir zeigen das auch in diesem Heft. 

Ich hoffe, dass Sie bereits ab dem nachsten DNL 
DERIVE- und Tl-Dateien von einer eigenen Homepa- 
ge abholen konnen. 


Dear DERIVE- and Tl-friends, 

Certainly you will be astonished at the photo- 
graph. No, I haven't fallen into any fountain of youth. 
May I introduce one of the youngest members of the 
DERIVE- and Tl community: Kimberly, our first 
granddaughter in the tender age of one hour. This 
picture will answer many questions (Hello, Barbel & 
friends). 

Here is now our 25th issue and 
on this occasion I’ve asked SWHH 
and SWHE for a contribution. They 
have both sent greeting addresses 
(page 37), many thanks. It is a pleas- 
ure for me to print them, but I’d like 
to share their appreciation with all of 
you. 

Sometimes I receive minor 
complaints: ... now and then the 
problems - and the accompanying 
solutions - seem to be from another 
star, very far away from my daily 
"Problemchens" with mathematics ... 
(D.BIum) 

Please consider that we fortu- 
nately have a big number of non-teachers in our 
group and that even the teachers should be allowed to 
reach for the stars now and then. I’d like to encour- 
age all of you who are feeling similar to report about 
their "Problemchens" (= small problems). No ques- 
tion, no letter, no contribution can be too simple or 
too unexciting. You all produce the DNL. And be- 
lieve me; I have the same maths "Problemchens" with 
the 15-19 years' students as you do. 

I have tried to make this DNL especially useful 
for the teachers' workaday routine. In the next DNL 
you will again find geometric items: mappings of 
objects in R 3 , hidden lines, fractals etc. There is a box 
full of Tl-questions with answers for you waiting to 
be opened. 

I think that in many cases DERIVE contribu- 
tions might inspire you for Tl-applications and vice 
versa. You will find examples in this issue. 

I hope that I will be able to offer the facility to 
download DERIVES and Tl-files from my school’s 
home page in the near future. 



Bis zum nachsten Mai 
Josef 



Until the next DNL 
Josef 



(*) If you are interested about the change of appearance of Kim between 1997 and now then 
have a look on page 48. There you can find a picture of Kim together with her 3 sisters and 3 
cousins who accomplished our unique “set” of seven grandchildren, Noor and Josef (proud 
grandparents). 
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Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
DNL. The more contributions you will 
send the more lively and richer in contents 
the DERIVE Newsletter will be. 


Austria 

Phone: 43-(0)660 31 36 365 
e-mail: nojo.boehm@pgv.at 

Preview: Contributions for the next issues 

3D-Geometry, Reichel, AUT 
Algebra at A-Level, Goldstein, UK 

Graphic Integration, Linear Programming, Various Projections. Bohm, AUT 

A Utility file for complex dynamic systems, Lechner, AUT 

Examples for Statistics, Roeloffs, NL 

Linear Mappings and Computer Graphics, Kiimmel , GER 

Solving Word problems (Textaufgaben) with DERIVE, Bohm, AUT 

Line Searching with DERIVE, Collie, UK 

About the "Cesaro Glove-Osculant", Halprin, AUS 

Hidden lines, Weller, GER 

Fractals and other Graphics, Koth, AUT 

Experimenting with GRAM-SCHMIDT, Schonefeld, USA 

Implicit Multivalue Bivariate Function 3D Plots, Biryukov, RUS 

The TI-92 Section, Waits a.o. 
and 

Setif, FRA; Vermeylen, BEL; Leinbach, USA; Halprin, AUS; Speck, NZL; 
Weth, GER; Wiesenbauer, AUT; Aue, GER; Propper, GER; Koller, AUT; 
Stahl, USA; Mitic, UK; Tortosa, ESP; Santonja, ESP; Wadsack, AUT; 
Schorn, GER, Chaffee, USA and 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: MagJosef Bohm 
Herstellung: Selbstverlag 


The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE User Group. It is pub- 
lished at least four times a year with con- 
tents of 40 pages minimum. The goals of 
the DNL are to enable the exchange of ex- 
periences made with DERIVE as well as to 
create a group to discuss the possibilities 
of new methodical and didactical manners 
in teaching mathematics. 

We include now a section dealing with the 
use of the TI-92 and we try to combine 
these modem technologies. 

Editor: Mag. Josef Bohm 
A-3042 Wurmla 
D'Lust 1 
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Problems: 

- Duffs awk can run out of memory with large files (more than a 10x10 grid, say, depending on how much 
low memory you have available.) 

- Classic DERIVE does not leave enough available memory to oad mawk through "Options Execute" (use 
Duffs awk or run ACRO.BAT outside DERIVE). DERIVE XM is OK. 

- Other awk interpreters may or may not have enough memory to run within DERIVE (gawk does not). The 
old large memory model compilation of mawk 1.0 (bmawk) does slightly better that the new version, but 
does not seem to be currently available from the net. 

Additional notes: 

- You can save a few keystrokes by faking ACROSPIN.EXE through the following C miniprogram 

#include<stdlib . h> 
void main ( ) 

{ system ( "aero . bat " ) ; } 

Compile as ACROSPIN.EXE and place in your DERIVE directory. Just use "aero" as the save file name. 

- the 3dv data file can be used for producing high quality hardcopy on a variety of devices. Look into 
3DVKIT1.ZIP, found whereever you found 3DV25.ZIP. 

- The memory problems could be solved if somebody with access to an awk compiler would post a compiled 
version of aero. awk. Or write it in C or whatever (but be careful to allow for multiple superimposed plots in 
the file.) 

Happy spinning! 

DNL: I followed Oscar Garcia 's hints and tips. And it worked. Fortunately some times ago R. Schorn 
had sent 3dv.exe, but I also couid find it on the net (together with KIT. ZIP and MAWK). I have contact 
to Oscar and it seems that he is trying making other DERIVE files of 3D-objects suitable for his 
3dviewer. The files produced by my ACD.EXE (DNL#24) have another format than the DERIVE- 
produced .ACD-files. Oscar sent also a patch to include automatic scaling. If you are interested in that , 
then please call or write, 111 send it to you. I add a screen shot of 3dviewer: 



(3dv is still working the DOS-environment, Josef) 


John Alexiou, USA j a72@prism.gatech.edu 

I am a Mechanical Engineering graduate student at Georgia Tech. We are using DERIVE for 
WINDOWS for Multibody Dynamics and I was wondering if there are other people out there doing 
the same. I was interested in exchanging some information and maybe .MTH-files. 

I have noticed that in many packages people use a LIM ( f (x) , x -A a) for a substitution, unfortu- 
nately this is SLOW and cannot handle vectors very well. 

I have a solution to this problem. You can use the ITERATE function to iterate once all the variables 
with their values. For example 
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#1: SUBEQ (eq, expr) := ITERATE (eq, LHS (expr) , RHS (expr) , 1) 


User 


# 2 : 


#3: 


SUBEQ 


SUBEQ 


- r x 


V (l 

- x-y) 

2 

t 

2 

LL X + 

y 



- r x 


V (l 

- x-y) 

2 

r 

2 

LL X + 

y 




t [x — If y — 0] 


= [ 1 / 1 ] 


f [x = a + b, y = a - b] 


User=Simp (User) 


User 


#4: 


a + b 2 2 

, V(- a + b + 1) 


Simp (#3) 


L 2 • (a + b ) J 

It applies to any vector/matrix size/shape, and to number/symbolic computations. In fact you can im- 
prove RK and EULER in ODEAPR. 
john 


Terence Etchells, Liverpool, UK t.a.etchells@livjm.ac.uk 

Hi ALL 

Eve done a fair bit of programming in my time using DERIVE , and to my knowledge there are no 
books on it at the moment. However, I have a book on Numerical Analysis published this month. 
There are many activities in this book that deal with programming DERIVE ; for example program- 
ming the Newton Raphson Method; Cobweb diagrams; Secant method; Lagrange interpolating poly- 
nomials; Simpson's rule; Richardson's method, Euler and Runge Kutta methods for ODEs etc. etc. 

Contact Philip Yorke at Philip@chartwel.demon.co.uk for further details. 

Incidently. 

I have been playing around with operators recently with DFW 4.03 and I discovered for myself some 
interesting mathematics. Try this: 

The operator D on f(x) we define as Df(x) = f(x+h) - f(x). (I would prefer to use capital delta for this 
but I doubt whether the character will transmit through e mail). 


#1: F(x) := 

A(p P x, h) := ( lim p) - p 
#2: x->x + h 

Simplify A(F(x), x, h): 

#3: A(F(x) p x, h) = F(x + h) - F(x) 

Now I want the operator A to act on F again i.e. A A 2 F(x): No Problem 
#4: A(A(F(x) , x p h)) = F(x + 2-h) - 2-F(x + h) + F(x) 

What about A A 3 F(x) or even A A 10 F(x)7 We can use the ITERATE command 
#5: ITERATE(A(s , x, h), s, F(x), 3) 

What do we get? 

#6: F(x + 3 -h) - 3-F(x + 2-h) + 3*F(x + h) - F(x) 

#7: ITERATE(A(s , x p h) p s p F(x) p 10) 

Are you surprised that we get 

#8: F(x + 10- h) - 10 • F(x + 9-h) + 45-F(x + 8-h) - 120-F(x + 7-h) + 210.F(x + 6-h) - 252-F(x + 

5 -h) + 210- F(x + 4*h) - 120-F(x + 3-h) + 45-F(x + 2-h) - 10-F(x + h) + F(x) 


As you can see, we can use now the A, Josef. 
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Mr Schmidt sent a wonderful MATHEMATICA - graphic of a snail house together with the 
functions to produce it. I tried with DERIVE and brought the snail into life using ACD and 
ACROSPIN. Many thanks, Mr. Schmidt for your idea. The snail is a new object in my collec- 
tion of animated graphs. (The 1997 snail can be found on page 36.) 


[Notation := Decimal, NotationDigits := 3] 


TI 

1 

' cf )2 

pi 10 :i , 

10 

r := 0.0972- 

3 +0.125 

,20 J 


zO := 0.0465 


az := IF 


.r, - 

l 20 J 


<t> - pi 10 


fa = 1, 1, IF 

r 

FLOOR 

' «i> ' 

t 

<l> 

i 

1.1. I.sll 

l 


, 2 > 


2 

JJ 


HAUSO, 6) := VECTOR(VECTOR([COS((|>-pilO)-r-(l. 1 + COS(0)), SIN(<|>-pi 10) • r • (1. 1 + COS(0)), 
l’-SIN(0) + ZO], (I), 0, 57), 0, 0, 2-77, pilO) 


ZACKO, a, z) := VECTOR(VECTOR( [COS((|> -pi 10) • r • (1. 1 + az-COS(z-pilO)), SIN(<|>-pil0).r-(l. 1 + 
az.COS(z.pilO)), r-az-SIN(z-pilO) + zO], 0, 57), a, 1, 1.2, 0.2) 


[HAUSO, 0), V ECTOR (ZACK(<|>, a, u), u, -2, 5)] 




The respective Graphs produced in the 3D Plot Window of DERIVE 6 
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Terence Etchells, Liverpool, UK 


t.a.etchells@livjm.ac.uk 


DERIVE: The Next Generation 

You may or may not be aware that I am studying for a PhD in Neural Networks. ( As an aside I am 
writing an article for the DNL on how I programmed Derive to train a neural network.) My experi- 
ences lead to the following request. 

As you know DERIVE is my number one program. If I were to have only one program on my desert 
island PC it would be DERIVE version 4.03. As a teacher of 16-19 mathematics it has always been 
the best and the Windows version is superb. 

However, now as a user at the University level I can see why some universities prefer other CAS such 
as MAPLE(MATLAB) and MATHEMATICA. Granted they are more difficult to learn and to use; 
and are not very intuitive. But, they have an important facility that DERIVE does not. That is, what I 
call, Vertical Programming, i.e. 

For while 

If.... Then 

End While 

as opposed to, what I call, horizontal progamming, 
vector(if(.... then ),r,...) 

I have attatched a MTH file that performs Neural network training. The effort required to program 
Derive to do this was considerable, I really enjoyed doing it but thats by the by, and most of the effort 
was in eliminating repetition of calculation through auxiliary functions. Something that I am sure, 
although I don't know, would be made easier with vertical programming. 

Another aspect which I find a little frustrating is the way DERIVE works as a Turing machine, i.e. all 
calculations are done before any ouptut is given. Which means that if I run a function that requires a 
long compute time (say 48 hours) and after 24 hours( which is required in my work) and wish to see 
where I am up to, I can't because interrupting loses all the information. Similary, having continous 
output (e.g. to look at convergence) as the function runs and producing dynamic plots of results that 
change with each new calculation is not possible. 

I am aware that a great deal of time and effort has been invested in the Windows version and as fan 
and power user of Derive I would ask (plead!) that a programming window be added in the next ver- 
sion. I think the programming window in the TI 92 is a good model on which to base a DERIVE pro- 
gramming window. I am sure, but again I don’t know, that all the programming functions such as 
FOR, WHILE, OUTPUT, IF, THEN, PLOT are already there in muLisp but the DERIVE user can't 
get at them or utilise them. 

I know that we users are always asking for additions and improvements and after such a mammoth 
effort as DfW more wishes must seem as a pain in the butt, but that is the nature of the beast. 

To end a little question, is a UNIX version of DERIVE a realistic possibility? (Please don’t swear! ! ! !). 

All my hopes! 
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Humberto M.Pereira Silva, Rio Tinto, Portugal 

Dear friend, 

. . . Reading newsletter #23 in which V. Hermans from Holland blamed for the lack of any 
incentive to pupils of High School to work with electronic machines in the classroom, will it be Texas, 
Casio, HP or else. In Portugal there are some new ideas two years ago when the use of electronic cal- 
culators were introduced in the last three terminal years (10, 11, 12) but unlikely were forbidden in the 
final examinations. I personally think that would be examination to test the skill of pupils in Calculus, 
because I think it will be more important in their future life than some subjects that has been taught in 
school benches. Im an electronic engineer who works and teaches the last form in High School and 
the first two years in University in Maths Analysis and Linear Algebra (matrices ...). 

At the other hand we have to congratulate for the excellent quality achieved with the new 
DERIVE for WINDOWS that I have been using for a month or so, it comes with some new fine fea- 
tures and is even easier to work quicker and powerful. 

Hoping to hear from you soon in meantime I wish you a nice and profitabel 1997. Yours 
H.M.P.S 


Soft Warehouse Hawaii swh@aloha.com 

Hello Power Users of DERIVE 

To keep you up to date, the following is a summary of the enhancements made in DERIVE for Win- 
dows version 4.03 

1. The quality of printed screen images of 2D and 3D plots has been improved by the elimination of 
spurious points. 

2. The date/time strings on printouts now use the format (dd/mm/yy or mm/dd/yy) appropriate for 
the user's own country. 

3. Comments in the DMO files are now handled like they were in the DOS version of DERIVE. 

4. A File Change Directory command is also included in the 2D and 3D-plot windows. 

5. In 3D-plot windows, it is no longer required that an expression be plotted before issuing the Write 
to Acrospin command. 

6. If ACROSPIN.EXE is in the YDfW directory (as it should be), then it will automatically be found 
by the Write To Acrospin command. 

7. When multiple expressions are highlighted and copied to the clipboard by pressing Ctrl+C, the 
highlighted rectangle is no longer corrupted. 

8. In 2D and 3D-plot windows, allow pi to be entered in upper or lower case. 

9. Allow quoted strings and underscores to be used in function and variable names given in the De- 
clare Function Definition, Declare Variable Value and Declare Variable Domain commands. 

10. If a vector is highlighted and the solve icon is clicked on, the Solve System instead of Solve Al- 
gebraically command is called (this was a serious deficiency in earlier versions of DfW). 
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The following enhancements are in both DfW and DfD version 4.03: 

11. A new built-in LOAD (file-name) function makes it possible to have DMO files automatically 
load utility files. 

12. The new function DIF_NUMERIC(y,x,xO,h,n) defined in NUMERIC.MTH numerically calcu- 
lates the nth derivative of y(x) at x = xO using a centered finite difference with a step of size h. 

13. The DSOLVE1 function in ODE1 .MTH is now better able to solve separable ODEs. 

14. Some problems with the DIRECTION_FIELD function defined in ODE_APPR.MTH have been 
corrected. 

15. Recognize that SUM((ax+b) A n,x,c,inf) can be expressed in terms of the Riemann Zeta function 
even if n is symbolic provided that n > 1 . 

Horst Scheppelmann, Hameln, Germany 

1st Ihnen ein Buch zu DERIVE bekannt, das sich iiberwiegend mit vektorieller analytischer Geo- 

metrie befafit, wie sie z.B. mit dem Lehrbuch von LAMBACHER-SCHWEIZER aus dem Klett- 
Verlag in der Sekundarstufe II (Klassen 11-13) unterrichtet wird. 

Mit freundlichen GriiBen und Wunschen fur das Neue Jahr. Ihr H.S. 

DA/L: The question is directed to our German friends, but maybe there is anybody among 
you who does know a book dealing with vectorial analytical geometry for Secondary Schools 
(forms 11-13)7 

Michael S Mullen, Austin, TX, USA msmullen@tenet.edu 

Why can DERIVE XM 3.01 not solve 2 A x = x A 2? The plot shows all three solutions including the 
obvious one x = 2. 

I tried to declare x as a complex variable and I tried to expand the equation prior to solving, but 

DERIVE just yields 2 A x - x A 2 = 0. 

Any assistance is appreciated. 

DNL: There were several answers, all of them with similar contents: 

One example: I had the same experience on the TI-92 until I changed the mode to Ap- 

proximate. 

x 2 

#1: S0LUTI0NS(2 - X = 0, X, -5, 10) DERIVE 6: 

#2 . |- 2 #2 is the “simplified” result and 

#3 the “approximated” one. 


[- 0 . 7666646959 ] 


Mike Hammet, Greenville, SC, USA Hammet_Mike/furman@furman.edu 

I have just received and started using DfW. I want to show you something that is good for comparing 
its speed with other versions of DERIVE and for checking the speed of your computer also. It is the 
following math problem. 

2 A 101 - 1 (which expands quickly gives) 

2535301200456458802992406410751 (which factored slowly) 

7432339208719 * 341117531003194129 
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Several years ago I ran this on a 286/8 MHz machine and it factored the number in 12 hours, 20 min- 
utes cpu time (That was with DERIVE 1.6). I wondered if it could even do the problem before I finally 
let it run overnight and returned the next day to find that it worked. Since then I have tried it on some 
much faster computers with newer versions of DERIVE. 

Here are the cpu times for some other computers and versions of DERIVE. 

DERIVE for WINDOWS: 486/50 - 22.3 minutes 

486/100- 12.1 minutes 

DERIVE 3: 486/50 - 47.0 minutes 

486/100 - 25.4 minutes 

DERIVE XM 486/50 - 24.0 minutes 

486/100 - 13.6 minutes 

Here DERIVE 3 and DERIVE XM were running under DOS only with WINDOWS off. By the way, 
running DERIVE 3 under WINDOWS slows it down by a factor of only about 1 or 2 percent but 
DERIVE XM is a disaster in WINDOWS, which slows it down by about a factor of five. It takes XM 
over two hours to factor the number of running under WINDOWS on a 486/50! 

So it seems DfW is even a little faster than XM. 

This is a fun problem to run on DERIVE. I would like to see the speed of a Pentium 200 on it. If you 
feel it is worth it please feel free to share it with other users. 

DNL : Some answers: 

M.S.Mullen: Your test was a great idea. 2 A 101— 1 really throws DERIVE into a tissy. After a little 
brute force I found that 2 A 93 - 1 is a fairly good test with a much shorter duration. 

On my 486DX100 the expand - factor sequence took 14 sec for DERIVE XM and 1 1 sec for DERIVE 

loaded low. On my old 386 SX16 XM took 63 sec 

Have fun 

M.S.Mullen has a nice question in his mail: 

" Who is General Failure , and why is he reading Drive C?" 

Scott Guth: Hello DERIVE enthusiasts! 

I thought Ed offer my results to FACTOR(2 A 101-1, Rational). 

Using DERIVE for WINDOWS: Pentium 90MHz/32MB RAM w/ WIN95: 13.6 minutes 

Using Maple V Release 3: Same computer as above; 4.02 minutes 

That sure says something about the Maple kernel! 

Dr.N.B.Backhouse: Dear Derivers! Let me give you some more timings of factorisations of 2 A 101— 1 : 

All done on a 133 Mhz Pentium (32 Megs of RAM, which may be or may not be relevant): 

DERIVE 3.00: 13.2 minutes; 

Macsyma 2.10: 9.0 minutes; 

Mathematica 3.0: 2.5 minutes 

Maple V5 Release 4: 1 .75 minutes 

Ubasic: 0.75 minutes, Didn't they do it well? 

Now can someone tell me if it is significant that both DERIVE and Macsyma are Lisp based, whereas 
Maple and Mathematica are written in C. (Read Al Rich s interesting answer on page 46, Josef) 
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Leon Magiera, Wroclaw, Poland 

Given: 

l 

f := 


magiera@rainbow.if.pwr.wroc.pl 




f 


2 A 

#1: 


2.E- I 

( VO ■ ro ^ 

2 k ■ r 


L 1 

l r J 

m j 



#2: R1 := j 


E 

— + 4 

V m 


#3: 

R2 := 4 4 


^ m 

#4: 

g(r) := J f dr 


E 

vv k j 
E 

LL k J 


2 2 » 
2 Vo ■ ro ■ m 


2 2 » 
2 Vo ■ ro ■ m 


DERIVE evaluates g(r):=INT(f,r) and g(R2)-g(Rl) but doesn’t simplify INT(f,r,Rl,R2). 
Best regards, Leon Magiera 


A. van der Meer, Twente, Netherlands 

Consider the following DERIVE session: 
l 


A.W.J.vanderMeer@math.utwente.nl 


#1: FOO := 


dx 


a + b- SIN GO 
Antiderivative calculated "by hand": 


#2: GOO := 


2 2 
4i a - b ) 


■ AT AN 


a- TAN 




+ b 


2 2 
4 (. a - b ) J 


d 1 

#3: — GOO = 

dx b-SINO) + a 

#4: G(2.tt| - G(0) = 0 

Conclusion: G is not the antiderivative on the whole real axis. In this case G has discontinuities on [0, 2*pi] 
For this example Derive has a trick to remove these discontinuities by adding a step function: 

Calculated by DERIVE: Simp(#l) 


2- AT AN 


# 5 : 


2-b-COSOO 


FOO := 


2-b-SINGO + G/(a - b) + 7 (a + b)) ^ 


7( a + b )-^(a - b) 


Which gives the correct results: 


2 -TT 


#6: F(2.jt) - F(0) = 


7(a + b) -7( a - b) 


2 -77 


# 7 : 


2 -7T 


a + b- SIN GO 


dx = 


V(a + b) -7C a - b) 


4 .jt 


#&: F(4.tt) - F(0) = 


7(a + b) -7Ca - b) 


(All screenshots are from DERIVE 6.10.) 


D 
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So, DERIVE is cleverer then you expected. In your example DERIVE has found that there is a possibility that 
certain combinations of values of the parameters may result in discontinuities in the antiderivative g(r) on the 
integration interval, but did not have a trick to remove these. I expect that the definite integral will be correctly 
calculated if you use realistic numerical values for the parameters. 

Klaus Fischer, Darmstadt, Germany 

Klaus Fischer suggests to include the most important and interesting WWW-sites for DERIVE and 
TI-92 users into the information page. Excellent idea , we will do that. Any recommendations are ap- 
preciated. See the Information page. 

Jos Verhoosel, Eindhoven, Netherland J.C.M. Verhoosel@pth.nl 

Dear DERIVE User, 

I teach mathematics for students who will become teacher in mathematics themselves. As part of their education 
they invest in their 2nd year 80 hours of programming (programming structures, recurrent programming, proce- 
dures, functions, subranges etc.) 

Recently I was asked if I was willing to replace this course with a new one: ’’Programming in DERIVE ’. I think 
this is a fantastic opportunity for them to extend their (and my) knowledge of DERIVE. 

Do you know if there exists such a book or course? To my knowledge there is none. 

Forced by circumstances, the course should start 3th february 1997 and concerns 2nd years' students in mathe- 
matics and I really would like to give the course. I already found some interesting examples in some utility files 
of DERIVE and in some Titibits - files from Johann Wiesenbauer. 

My question for you: 

Can you give me some more programming problems in DERIVE. The problems should cover most of the pro- 
gramming facilities of DERIVE Also, the problems should be simple to grasp for my students and should be 
practical. (That means that it should be utterly clear that such a utility or subprogram can be of great use during 
the rest of course) 

Can you give me some suggestions? 

(I intend to share my experiences in the International DERIVE Journal if interesting) 

Thanks in advance, Yours sincerely 
Drs. Jos C.M.Verhoosel 
Fontys PTH, Dept of Mathematics 
Eindhoven, The Netherlands. 

DNL: I could imagine that the "TURTLE GRAPHICS" article in this issue could contribute for your 
course and Terence's nwe book, of course. And now, once more Terence Etchells: 


Terence Etchells, Liverpool, UK t.a.etchells@livj m.ac.uk 

Hi Josef, 

I have written a little set of functions that will produce an interpolating polynomial to approximate a function at 
any given points on that function. So for example, you may wish to find a polynomial that passes through the 
points (0,sin(0)), (pi/2, sin(pi/2)) and (pi,sin(pi)). So our function is sin(x) and the x values of the desired points 
are [0, pi/2, pi]. As we have 3 points the polynomial will be a quadratic. You will find attatched a MTH file 
INTERPOL.MTH that will do this automatically. In fact it will work for up to 20 points (i.e. a polynomial up to 
order 19), this can be increased by increasing the available variables in the def vars (and the line above that 
resets all variables we are to use to empty) definition. 

Have you found an easier way (is there? A1 or Dave) of producing a vector of undefined variables, without hav- 
ing to type them all in. Something like VARIABLES(30,a) which would simplify to 

[a,b,c,d, ,x,y,z,aa,bb,cc,dd]. 
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You will also notice that on line #9 the LIM function variables vector needs an EXPAND command, it doesn't 
work without it. Any ideas why? 


If nobody has come up with a function that will do this before ( don't remember seeing one in any previous 
DNL) you might like to use it a forth coming DNL. 


#1: def_vars := [a, b, c, d, e, f, g, h, i, j, k, 1, m, n, o, p, q, r, s, t] 

#2: [FUNCTIONGO :=, rvalues :=] 

#3: PO L YNOM I AL_AUX (dim) := VECTOR^def_vars -x , r, 1, dimj ■ VECT0R(1 , r, dim) 

#4: POLYNOMIAL (x) := POLYNOMIAL_AUX(DIMENSION(xvalues)) 

#5: POLYNOMIAL_AUX(x , dim) := VECTOR^def_vars -x , r, 1, dimj ■ VECT0R(1 , r, dim) 

#6: POLYNOMIAL (x) := POLYNOMIAL_AUX(x , DIMENSION(xvalues)) 

vars_vect := VECTOR (def_vars , r, 1, DIMENSION(xvalues)) 

#7: r 

co_eff := (SOLUTIONS (VECTOR (FUNCTION (x values ) = POLYNOMIAL (xvalues ), r, 1, DIMENSION(xvalues)) , vars_vect)) 
#S: r r 

interpolate := lim POLYNOMIAL (x) 

#9 : EXPAND(vars_vect)-»co_eff 

Example: Interpolate cos(2x) + sin(3x) by a polynomial atx = 0,1 .S.SAS.SJ.S 


#10: 

FUNCTION(x) := C0S(2-x) + 

SIN(3-x) 

#11: 

xvalues := [0, 1, 2, 3, 4, 

5, 6, 7, 

#12: 

interpolate 


#13: 

TABLE(FUNCTION(x), x, xvalues) 

#14: 

TABLE(interpolate , x, 0, 

9, 0.005) 



DNL: Plot #12 and obtain the interpolating polynomial. #13 results in the nine given points. #14 would 
plot the interpolating polynomial in form of a thick line (point size medium or large and points con- 
nected). 
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AN IMPLEMENTATION OF "TURTLE GRAPHICS" 

IN DERIVE 3 

Josef Lechner 1 , Eugenio Roanes-Lozano 2 , Eugenio Roanes-Macias 3 , 

Johann Wiesenbauer 

^undesgymnasium Amstetten (Austria) 

2,3 Dept. Algebra, Univ. Complutense de Madrid (Spain) 

4 Institut fiir Algebra, Tech. Univ. Wien (Austria) 


Introduction 

If you have worked previously with any LOGO dialect you will know about the nice possibilities of 
the Turtle Graphics. We had already developed implementations for Turbo-Pascal [Ro4] and Maple V 
([Rol], [Ro2], [Ro3]). Now we have developed another one for DERIVE 3 [1] , that is presented below. 


Two Remarks 

To follow the article, a certain knowledge of the possibilities of Derive and Turtle Graphics is sup- 
posed. [Ab] is a very nice introduction to the later, and [A-dS] is an impressive collection of ideas for 
its use. 

Note that in some early versions of Derive 3 evaluation is made in such a way that this code does not 
work. The code has been checked in versions 3.13 (Classic and XM) and Derive for Windows Beta v. 
0 . 1 . 


Getting Started 


We should begin by setting the graphics in Connected Mode and the Color in the non Auto Mode. 
Then the Turtle file turtlesut.mth should be loaded normally (or as a utility file in the background). 

Example 0: For instance, to draw the square of vertices (0,0), (0,1), (1,1), (1,0) it is enough to Plot 
[ FD(1) , RT(90) , FD(1) , RT(90) , FD(1) , RT(90) , FD(1) , RT(90)] 



1.5 

+ 


0.5 


- 1.5 -1 - 0.5 

0.5 1 

- 0.5 

-1 

1.5 2 



m Works now also for DERIVE 6, Josef 
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Or, even more simple, to make the assignation 
to_rep:=[FD(l) ,RT(90)] 

(and to Simplify) and Plot 
REPEAT(4) 

Note that the syntax adopted for the REPEAT command is not the usual in Logo. 

Observe that, unfortunately, unlike in Logo, there is no visible cursor for the turtle in the graphics 
screen. 

About the implementation (the content of the turtlesut.mth file) 

The first two lines of the tu rtl es_ut . mth file are 

Angle:=Degree 

Preci si on :=Approximate 

(in order to increase the speed; nevertheless this can -should- be omitted in some cases -see [Ro3]-). 
The main problems we had to produce this implementation were related to the difficulties to change, 
from a programme, assignations made to variables. The following auxiliary functions are to be used 

NEWX(u) :=(xcor:=u) 

NEWY(u) :=(ycor:=u) 

NEWD(u) :=(headi ng :=u) 

The turtle is, at any moment at point (xcor, ycor) and heading towards ’’heading” (counting clockwise 
from the usual halfaxis y+). The user can check the value of these variables by simplifying them. 
Advancing a length 1_ can be implemented easily: 

FD(l_):=[xcor , ycor ; NEWX(xcor+l_C0S(90°-headi ng)) , NEWY(ycor+l_SIN(90°-headi ng)) ] 

and moving backwards the same distance can be obtained from the previous function 

BK ( 1 ) :=FD (— 1 ) 

When the turtle turns, only the direction (heading) is altered. The new value is reduced modulo 360 in 
order not to get high values for the ’’heading” variable after many rotations 

RT(a_):=[xcor,NEWY(ycor+ 0 • NEWD(MOD(headi ng+ a_°,360°)))] 

LT(a_) :=RT(-a_) 

The last of the typical turtle-commands is ’’home” (return to the centre of the screen facing upwards): 
home := [NEWX(O), NEWY(O) , NEWD(O)] 

But there are some more standard turtle-commands mixing turtle ideas and Cartesian coordinates. To 
go to a certain position of known coordinates can be done using 


SETPOS (x_ , y_):=[ xco r , y co r ; xco r:=x_ , y co r:=y_] 0 






D-N-L#25 


L&R&R&W: TURTLE GRAPHICS 


P 17 


The actual position is given by 
pos := [xcor, ycor] 


And to change only one of the coordinates of the current position can be done using 

SETX(x_) :=[[xcor,ycor] , [xcor:=x_,ycor]] 

SETY(y_) :=[ [xcor , ycor] , [xcor , ycor : =y_] ] 


It is also possible to force the turtle to look in direction a_ with 

SETH(a_):=[xcor ,NEWY(ycor + 0-headi ng:=a_) ] 

or to look towards a certain point (x, y) with 

SETHTOWARDS(x_,y_):=[xcor , NEWY(ycor+0-headi ng:=ATAN(x_-xcor ,y_-ycor)) ] 

Note that, to clear the screen, there is no ClearScreen command. The Delete/All in Derive's Plot menu 
should be used instead. 

Another standard feature of Logo is the REPEAT command. We couldn't maintain the standard syntax 
but the idea is similar. 

REPEAT(k_) := IF(k_ = 1, to_rep, APPEND(REPEAT(k_ - 1), to_rep)) 

I tried successfully: 

REPEAT(k_) :=VECTOR(to_rep, i, k_) 

Finally, the variables are initialised: 

xcor :=0 
ycor :=0 
headi ng :=0 


About REPEAT 

What has to be repeated has to be stored in the to_rep variable (list) and the only argument of REPEAT 
is how many times has to be repeated what is stored in that list. Therefore, to execute repeatedly a list 
of commands just type 

to_rep:= [ List of commands ] 

REPEAT( number) 

and Plot. 

Observe that, the way it is defined, a REPEAT can not be nested inside another REPEAT. Defining a 
function REPEAT2 depending on a variable to_rep2 does not always work. 
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Some Simple Examples 


Remember: 

i) To Simplify “home” and to choose Delete/All in the Plot menu in between executing the exam- 
ples. 

ii) (Simplify and) Plot each new expression (which is not an assignation). 


Example 1 : Draw an L. 

#2: home 

#3: [xcor := 0, ycor := 0, heading := 0] 

#4: [LT (90) , FD(2) , RT(90), FD(3)] 


Example 2: Draw a flag and return to the origin using only the Cartesian commands. 


#5: [xcor := 0, ycor := 0, heading := 0] 

#6: [SETY(2), SETX(2) , SETY(l) , SETXC0) , SETY(O)] 


Example 3: Draw a regular heptagon 


#7: [xcor := 0, ycor := 0, heading := 0] 


#8: to_rep 


FD(2) , RT 



#9: REPEAT (7) 



Example 4: Draw regular polygons of 5,7 and 9 sides of length 1 


#10: 

[xcor := 0 

#11: 

n_ := 

#12: 

to_rep := 

#13: 

n_ := 5 

#14: 

REPEAT (5) 

#15: 

n_ := 7 

#16: 

REPEAT (7) 

#17: 

n_ := 9 


ycor := 0, heading := 0] 


FD(1) , RT 



3 



#1S: REPEAT (9) 
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Example 5: Draw the 36 wires of a chart- wheel 


#19: [xcor := 0, ycor := 0, heading := 0] 

#20: to_rep := [FD(1), BK(1) , RT(10)] 

#21: REPEAT (36) 



PenUp and PenDown 

There is no PenUp/PenDown command in our implementation. Instead, the user has to play with only 
Simplifying (that only changes the position of the turtle) or Simplify and Plotting (that changes the 
position of the turtle and draws the correspondent segments). 

Example 6: Draw two segments, the first one with endpoints (0,1) and (0,3) and the second one with 
endpoints (1,3) and (3,3). 


#2&: [xcor := 0, ycor := 0, heading := 0] 

#29: FD(1) 


#30: 


0 0 
_ xcor := 0 ycor := 1 _ 


#31: FD(2) 


#32: [RT(90) , FD(1)] 


#33: 


[0 , ycor := 3] , 


#34: FD(2) 


0 3 1" 

xcor := 1 ycor := 3 


#29 and #32 must be simplified only; 
#3 1 and #34 must be plotted. 


3 

2 

1 


-1 


1 


A Recursive Example 

Example 7: Function KOCH below constructs the side of a koch star with length “len” and depth “n”. 


K0CH(len, n):=IF(n > 0, [K0CH(len/3, n - 1) , LT(60) , K0CH(len/3, n - 1) , 

RT (120) , K0CH(1 en/3 , n - 1) , LT(60) , K0CH(len/3, n - 1)], FD(len)) 


#36: [xcor := 0, ycor 0, heading 0] 

#37: RT (90) 

#38: [0, ycor := 0] 

#39: K0CH(1, 1) 

Simplify #37 (in order to have the x-axis as base) 
and then plot #39. 
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depth 2 and 4: 

#40: [xcor := 0, ycor := 0, heading := 0] 

#41: RT (90) 

#42: [0, ycor := 0] 

#4E : KOCH (1, 2) 



#44: [xcor := 0, ycor 0, heading 0] 

#45: RT (90) 

#46: [0, ycor := 0] 

#47: KOCH (1, 4) 



KOCH can be used to draw stars (= snow flakes): 


#4S: 

[xcor := 0, ycor := 

0, 

heading := 

#49: 

to_rep := [K0CH(1, 

4), 

RT (120) ] 

#50: 

REPEAT (3) 




In DERIVE 6 we can write a short pro- 
gram, Josef: 


kochstar(len , n, star, n_) := 

Prog 

star := to_rep 
n_ ;= 1 

#51: Loop 

If n_ = B 

RETURN star 

star := APPEND(star, to_rep) 
n_ : + 1 

#52: [xcor := 0, ycor := 0, heading := 0] 

#53: SETPOSCL. D 

#54: to_rep := [K0CH(1, 3), RT(120)] 

#55: kochstar(2, 3) 
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This is the original DOS-DERIVE screen from 1997: 



Conclusions 

Many geometric designs and drawings can be produced in a more convenient way using Turtle 
Graphics instead of Cartesian coordinates. Therefore, we think that this can be another useful tool for 
DERIVE 'rs. 
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Now is 2011 and we have DERIVE 6.10 which is much more powerful than the DOS-Version 
from 1997. 

I remember my LINDENMAYER-Systems contribution from DNL#51 and DNL#52. This has a 
tight connection to Turtle-Graphics. 

I use turtles_ut.mth for reproducing some examples: 
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hilb(len, n) := IF(n > 0, [ hi 1 b (1 en/3 , n - 1), LT(90), hilb(len/3, n - 1), 
RT(90), hilb(len/3, n - 1) , RT(90), hilb(len/3, n - 1) , LT(90), hilb(len/3, 
n - 1)], FD(len)) 


#53: 

0 

1 1 

1 

0 
u 

X 

1 1 

ycor := 0, heading := 0] 

#59: 

RTC90) 


#60: 

[0, ycor := 

0] 

#61: 

hi lb(2, 1) 




■0.5 



-0.5 

0.5 1 1.5 


#62: [xcor := 0, ycor := 0, heading := 0] 

#63: [0, ycor := 0] 

#64: hi lb(2, 4) 



pat2(len, n) := IF(n > 0, [pat2(len/2, n - 1) , LT(90), pat2(len/3, n - 1) , 
RTC180), pat2(len/3, n - 1) , LT(90) , pat2(len/2, n - 1)], FD(len)) 


#67: [xcor := 0, ycor := 0, heading := 0] 

#63: pat2C2, 1) 


2 

1.5 

1 

0.5 


-1 -0.5 



#69: [xcor := 0, ycor := 0, heading := 0] 

#70: pat2C2, 5) 
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#71: [xcor := 0, ycor := 0, heading := 0] 

#72: RT(90) = [0, ycor 0] 

#73: VECT0R([pat2C2, 5), LT(90)], i, 4) 


twig(len, n) := IF(n > 0, [twig(len/3, n - 1) , LT(45), twigClen/3, n - 1) , 
BKClen/3), RT(45) , twigClen/3, n - 1), RT(45) , twig(len/3, n - 1), BK(len/3), 
RTC315), twigClen/3, n - 1)], FD(len)) 


PTTmTmMTm 

I -Fh--4+- I 


4 
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#75: [xcor := 0, ycor := 0, heading := 0] 

#76: twig(3, 1) 


#77: [xcor := 0, ycor :- 0, heading := 0] 

#78: twig (3 , 5} 


More References: 



Lindenmayer Systems: DNL#51 andDNL#52 

Snowflake a.o.: DNL#33 and DNL#34 (Contribution from Matia Koth) 

LOGO in DERIVE: DNL#38 (Josef Lechner) 
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twig2(len, n) := IF(n > 0, [twig(len/3, n - 1) , LT(40), twig(len/3, n - 1) , 
BKClen/3), RT(45) , twig(len/3, n - 1) , RT(40) , twig(len/3, n - 1), BK(len/3), 
RT(310), twig(len/3, n - 1)], FD(len)) 


#30: [xcor := 0, ycor := 0, heading := 0] 

#31: twig2(3 , 5) 



twig3(len, n) := IF(n > 0, [twig(len/3, n - 1) , LT(30), twig(len/3, n - 1) , 
BKClen/3), RT(45) , twig(len/3, n - 1) , RT(30) , twig(len/3, n - 1) , BK(len/3), 
RT(250), twig(len/3, n - 1)], FD(len)) 
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#33: [xcor := 0, ycor := 0, heading := 0] 

#34: twig3(3 s 5) 





(Two more examples on page 28) 
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In an earlier Letter of the Editor I mentioned that the DERIVIAN Turtle is one result of the 
many fruitful talks at the DERIVE Conference 1996 in Bonn. (The most important are the 
breaks.) After this conference began a very busy exchange of ideas via e-mail. Following the 
TURTLE. MTH file you might ask yourself: 'what is that, there are two := assigments in one 
expression?. So did I. 

Fortunately Josef Lechner is very conscientious and he collects many things concerning 
Computer Algebra in general and concerning DERIVE in special. When I asked Josef to give 
an additional explanation for this strange expression he took two sheets of paper from his 
huge desk and asked me: "Will that be enough? I summarized my e-mails with Johann 
Wiesenbauer. " Many thanks, Josef. I think that will do it. So find here instead of his TITBITS, 
J.Wiesenbauers "double assignation". 


Comment of 201 1 : The following problems are from 1997. Later versions of DERIVE don't 
behave the same. I recommend experimenting with a := expression(a) applying the various 
possibilities of “Simplify” from the menu- and the edit bar as well !! 


Is it able to implement TURTLE graphics in DERIVE ? 

Not completely! 

The reason is that DERIVE does not support self assignments. In DERIVE we meet a strict functional 
programming style. Especially users, which are not accustomed with that way or programming might 
be frustrated by the result of the assignment 

a := expression (a) 

One example: The user enters 


#1: 

a : = 3 

User 

#2: 

a := a + 7 

User 

and expects a = 10, but he obtains a = a + 7, and simplifying again and again: 


#3: 

a + 7 

Simp (#2) 

#4: 

a + 14 

Simp (#3) 

#5: 

a + 21 

Simp (#4) 

So, what can be done? Let's try with the assignment NEWX(u) : = x := u with u being an expression 
in x. That works because expressions in function parameters are evaluated before passing. That is an 
(interesting) alternative to the ITERATES command. 

#7: 

NEWA(u) := (a := u) 

User 

#8: 

a : = 1 

User 

#9: 

NEWa[L— . j^ a + — -]] 

User 

#10 

1.5 

Simp (#9) 

#11 

a = 1.5 

User=Simp (User) 

#12 

1 .41666 

Simp (#9) 

#13 

a = 1.41666 

User=Simp (User) 

#14 

1.41421 

Simp (#9) 

#15 

a = 1.41421 

User=Simp (User) 
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Obviously the variable has changed its value! ! 

Let's try once more. Set back the value for a. Watch the annotations, please. 

#17 : a : = 1 User 

#18: 1.5 S imp ( # 9 ) 

#19: 1.41666 Simp(#9) 

#20: 1.41421 Simp(#9) 

#21: 1.41421 Simp (#9) 

#23: [NEWX(u) := (x := u) , NEWY(u) := (y := u) , NEWA(u) := (a := u) ] User 

#24: FWD ( 1 , w) := [NEWX(x + l-COS((a + w)-°)), NEWY (y + l-SIN((a + w ) - ° ) ) 

+ 0 • NEWA ( a + w) ] 

It seems to work: 

#25: home := [x := 0, y := (a := 0)] User 

#26: quadrat := [home, FWD(2, 0), FWD(2, 90), FWD(2, 90), FWD(2, 90)] User 


#27: 2 2 


S imp (#26) 


The big disillusion: 


#28: home = [0, 0] 


User=Simp (User) 


#29: VECTOR ( FWD ( 2 , 0), i, 1, 4) 


User=Simp (User) 


#30: InputMode := Word 


#31: ITERATES (FWD (2, 0), turtlepos, home, 4) 


User=Simp (User) 


2 0 
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What is possible in a list, does not work in a VECTOR! ! 

#33: a := 1 User 


r 

r 1 

r 2 i 



r 1 

r 2 i 



r 1 r 

NEWA 




, NEWA 




, NEWA 



a i - 



a t 



- 

- 2 

a - 



- 2 

a - 



- 2 


#35: [1.5, 1.41666, 1.41421] 

#36: a : = 1 


#37: VECTOR 

NEWA 

r 1 

r 2 i 

a + 

, i, 1, 3 



- 2 

a - 

- 



User 


Simp ( #34 ) 
User 


User=Simp (User) 


Conclusion: To implement a Turtle graphic in DERIVE it would be necessary for the program to sup- 
port a self assigment of the variables and/or the the support of composition of functions would be im- 
proved. 


So I will close: And yet it works. I tried to translate Josef (L) 's comments. I hope that 
he will recognize the sense of his comments in my words. Josef (B). 

J. Wiesenbauer has promised to check the DERIVE utility files for improvements 
using the new built in DfWand DfD4 functions and capabilities. Do you have ideas? 
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Tilgung fremderregter Schwingungen durch 
abgestimmte Ankopplung 

Cancellation of Separately Excited Oscillations 
by a Tuned Connection 

Leo H. Klingen, Bonn, GER 

Eine rotierende Maschine mit exzentrischen beweglichen Teilen steht auf einem Fundament, das beim 
Lauf der Maschine moglichst in Ruhe bleiben soli. Nicht immer stellt eine moglichst groBe Masse des 
Fundaments eine optimale Losung dar. Bei schwingungsfahigem Fundament mi haben wir den wohl- 
bekannten Fall der Federschwingung (Federkonstante Ci) mit periodischer Fremderregung x e , die hier 
beim prinzipiellen Ansatz von Figur 1 am Aufhangepunkt der Feder angreift. 

A rotating machine with eccentric moveable parts is 
placed on a basement which should keep as motionless 
as possible when the machine is running. Sometimes the 
optimal solution is not to have the mass of the foundation 
as big as possible. With an oscil Hating (vibrating) base- 
ment m-i we recognize the well known case of a spring 
oscillation (spring constant c-i ) with periodical separate 
excitation which acts in the spring's suspension point 
(figure 1). 

Due to Hooke's Law and Newton's Dynamic Fundamen- 
tal Law we obtain the motion equation when neglecting 

damping. 



Figur 1 

Mit den Daten Ci = 1, mi = 2 und x e = 4 cos (V3 t) und der realen Federdehnung (xi - x e ) ergibt sich 
die Bewegungsgleichung nach dem Hookeschen Gesetz und Newtons dynamischen Grundgesetz, 
wenn Dampfung vernachlassigt wird, zu 

2*i = -(*! -* ) 

oder/ or 

1 t . x +0.5x = 2cos(V3n 

Man kann sie exakt losen mit V / 

#3: DSOLVE2_IV(0 , 0.5, 2.COS(73-t), t, 0, 0, 0) 


#4: 










p 30 


Leo Klingen: Schwingungen / Oscillations 


D-N-L#25 


Oder man wendet nach Substitution ,Vj = .V3 die Runge-Kutta-Methode an: 

#4: (RK([x3, - 0. 5-xl + 2-COS(V3-t)] , [t, xl, x3], [0, 0, 0], 0.3, 50));r[l, 2] 

#5: 4 ■ COSC/3 • t) 


Das Ergebnis zeigt Figur 2 mit einer erheblichen Schwingungsantwort des Fundaments. 

We can find an exact solution using DSOLVE2Q or an approximative solution using the 
Runge-Kutta-Method. The result in figure 2 shows a serious vibration response of the base- 
ment. 

We can achieve an interesting extinction of the resulting oscillation by connecting an addi- 
tional mass m 2 using a second spring (spring constant c 2 ). See figure 3. 



Figur 2 


Eine interessante Ausloschung der resultierenden Schwingung kann man durch Ankopplung einer 
weiteren Masse m 2 iiber eine zweite Feder mit der Federkonstanten c 2 erreichen. Vgl. die Prinzipskizze 


Figur 3. 



Mit den zusatzlichen Daten m 2 = 1, c 2 = 3, die so gewahlt sind, dass 
sich als Eigenfrequenz die Frequenz des Storers ergibt, erhalt man das 
System (unter Beriicksichtigung richtiger Vorzeichen nach dem Prin- 
zip actio = reactio) 

The additional data (m 2 = 1, c 2 =3) are chosen in such a way 
that the proper frequency results from the frequency of the dis- 
turber. (principle of actio = reactio). 

Then we use again Runge-Kutta. 


oder 


x\ = - 0.5(xj - x e ) - 1.5(xj - X 2 ) 
x 2 = -3(x 2 -xj) 


Xj + 2x\ - 1.5x2 


2 cos(V3 t ) 


x 2 - 3xj + 3x 2 = 0 


Figur 3 
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Mit den Substitutionen ij = X 3 und %2 = *4 kann man Runge-Kutta ansetzen. 

#6 : (RK([x3 , x4, - 2-xl + 1.5-x2 + 2-COSG/3-t) , 3-xl - 3-x2], 

[t, xl, x2 , x3 , x4] , [0, 0, a, 0, 0], 0.3, 50))u[l, 2] 

First of all we have to find the initial values for a (that is the constant amplitude of the con- 
nected oscillation). We plot for two values: a 1 = -2/3 and a 2 = -4/3 (negative, because oppo- 
site to the spurious oscillation). Using a 1 we recognize a remaining vibration whereas using 
a 2 the effect of the additional oscillator causes a complete extinction of the basement's vi- 
brations. 

Zuvor muss der Anfangswert a (das ist bei vemachlassigter Dampfung die konstante Amplitude der 
angekoppelten Schwingung) bestimmt werden. Wir plotten fur zwei Werte ai = -2/3 und a 2 = -4/3 
(negativ wegen der Gegenphasigkeit zur Storschwingung). Beim Wert ai bleibt, wie die Grafik verrat, 
eine Restschwingung des Fundaments iibrig. Beim Wert a 2 ist die Gegenwirkung des Zusatzschwin- 
gers so groB, dass sich eine vollstandige Ausloschung der fremderregten Schwingung des Maschinen- 
fundaments ergibt! 


•x2 + 2-COS(V3-t), 3-xl - 3-x2], [t, xl, x2, x3, x4] , 0, 0, 

Utl, 2] 


•x2 + 2-COS(73-t), 3-xl- 3 • x2 ] , [t, xl, x2, x3, x4] , 0, 0, 

III 1 , 2 ] 

•x2 + 2-COS(73-t), 3-xl- 3-x2] , [t, xl, x2 , x3, x4] , 0, 0, 

U[l. 3] 

Natiirlich ist die Bestimmung von a nicht nur durch Plotexperimente, sondem auch rechnerisch leicht 
moglich. Man unterstellt naheliegend x, = X l cos ( v"" 3 1) bzw. x 2 = X 2 cos (V3 1) als partikulare Losung 
und berechnet fur A', - 0 fur alle t 

X 2 = a 2 = -4/3. 
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Figur 4 


Fur die notwendigen zweiten Ableitungen und algebraischen Umformungen hilft DERIVE , Rechen- 
fehler zu vermeiden. Ubrigens ist auch ein Realexperiment ohne sonderlichen Aufwand fur Liebhaber 
oder im Rahmen eines Schiilerprojektes ohne weiteres moglich. 

Ein elektrisches Analogon besteht in einem Antennen-Schwingkreis mit (mehreren) von Sendern 
fremderregten Schwingungen, darunter vielleicht eine besonders starke eines Lokalsenders. Wenn 
diese unerwiinscht ist, kann man sie "heraussaugen" durch Ankopplung eines abgestimmten zweiten 
Schwingungskreises mit eventuell verstarkter Frequenz des Lokalsenders, sodass das Frequenzge- 
misch des Restes ungestort in der iiblichen Weise weiterverarbeitet werden kann. Man kann versu- 
chen, auch diese Situation zu simulieren. 

An electronic analogon is an antenna oscillation circuit with (some) undesirable frequen- 
cy(ies) - local station. You could "suck it out" by connecting a second circuit with possibly 
amplified frequency. You could try to simulate this situation. 
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Probability Distributions 
Proof and Computations (2) 

Peter Mitic, Medstead, UK 

MEANS AND VARIANCES OF CONTINUOUS RANDOM VARIABLES 

Similar problems are encountered with continuous distributions, for which the moment gen- 
erating function is M(t) = J e xt f(x)dx where f(x) is a probability density function. When cal- 

culating the mean and variance of some continuous random variables it appears that 
DERIVE is more successful in evaluating integrals directly rather than via a moment generat- 
ing function. The DERIVE session below shows that integrals which do not involve the pa- 
rameter t can be evaluated easily. Integrals which contain this parameter cannot always be 
evaluated, even if the range of the parameter is restricted such that the integral is convergent 
(t < 1/2 in the expressions # 5 -.# 7 below). In particular, it is easy to check that the area 
under the graph of a probability density function is 1 (expressions # 2 to # 4 below). The 
probability density function used (F(x) in expression # 2) is that of a x 2 (5) random variable. 

- x/2 3/2 

</2 ■ e ■ x 

#1: F(x) := 

6 ■ J7T 

#2: Check that the area under the graph is 1 



# 8 : 


Failed! Calculate mean and 
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#9: mean := J X'F(x) dx 

0 

#10: mean := 5 


#11: variance := 


■ F(x) dx - mean 


#12: variance := 10 

The situation with respect to parameters is not at all consistent. In the case of a Normal(0,1) 

e 2 

random variable, with probability density function N(x) = , , DERIVE can immediately 

V 2;r 

simplify the expression for the moment generating function, 
M(t) :=int (N (x) *#e A (x*t) ,x, -inf , inf ) to obtain the result #e A (t A 2/2 ) with no 
apparent need to restrict the range of t. This result is required for a proof of the Central Limit 
Theorem. (See the next DNL, Josef). 

We have already pointed out (Mitic and Thomas 1994) that problems arise when the com- 
puter fails to calculate a result, and that consequently, problems presented as student exer- 
cises must be chosen so that they are 'do-able'. This is the case here. In order to teach mo- 
ment generating functions, one has to choose suitable probability density functions. We have 
found that probability density functions of the form g(x) x n , x> 0 with n a positive integer work 
provided that n is explicit. General results may be conjectured by constructing a VECTOR of 
particular examples. In the context of the class of probability density functions 


f(x,n) = 


x eZ + , the expression 


VECTOR (INT (x A n #e A (-x)/n!, x, 0, inf), n, 0, 10) 

can be used to check that the areas under the graphs of the first 1 1 members of this class 
are all 1. 


We would prefer to see the general result (obtained using MATHEMATICA 2.2, which does 
not simplify the result to 1) 


dx = 


r(l + n) 


into which we can substitute n = 0, 1 , 2, ... 


The need to consider particular examples in this way is regrettable, but it is unavoidable with 
DERIVE. 
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CALCULATING PROBABILITIES 

Calculating probabilities presents less of a problem technically, and also reinforces concepts 
in a way that a statistical package cannot. 

A frequent task with discrete distributions is to use the probability generating function, G(r), to 
calculate probabilities by isolating a coefficient or coefficients of t in a series expansion of 
Git)- This often causes computational problems which overshadow the point of the exercise. 
We illustrate (#14 - #16 below) how this task is very simple using DERIVE. The geometric 
distribution Geom(1/6) has the probability generating function given by #14. The coefficient 
of t r in G(t) is obtained by computing 

G(t) 

No formal expansion of G(t) is needed, so the student can concentrate on the precise mean- 
ing of terms in a probability generating function. If the substitution t = 0 is programmed, as 
below, it is necessary to use the limit construct shown, despite its technical incorrectness. 

t 

#14: G(t) := 

6 - 5-t 

#15: For obtaining P(X = 3): 

fd >3 
- G(t) 

#16 : UltJ 25 

lim = 

t-*0 3 ! 216 

In the case of continuous random variables it is best to use a probability density function di- 
rectly. The use of a direct integration in this way stresses the concept that a probability is 
represented by the area under the graph of a probability density function, whether or not the 
function concerned is integrable analytically or not. 

Expressions #18 -#20 below illustrate the calculation of P(0 <X<2) where X ~ j 2 (5) . 

-1 

2 e ■C3->-e-ERFCl) - 10) 

#13: J FGO dx = 

0 3 

#19: simplified above and the right side approximated below: 

2 

#20: J F(x) dx = 0.1503549639 

0 
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CONCLUSION 

The point of using a computer algebra system should really be that any problem can be tack- 
led. This context has highlighted a class of integrals and series which, unfortunately, cannot 
be tackled successfully with DERIVE. However, given certain general results, DERIVE can 
be used to compute means, variances, moments, probabilities etc without algebraic difficul- 
ties. Furthermore, all of these can involve some symbolic parameters, so that functional rela- 
tions can be explored. This is not possible with a package which can only manipulate num- 
bers. The shortcomings of DERIVE can be turned to advantage if the student is constrained 
to study relevant techniques and concepts, and hence to do some parts of the computation 
by hand. 
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A Toast Celebrating the 25th Issue of the DERIVE® Newsletter 

Josef and Noor, we in the DERIVE community are all so very grateful for your inspi- 
ration to initiate the independent DERIVE User Group and its bulletin; and for your talent and 
dedication in making it work so well. 

Many of us are surely thinking “Can it really be 25 issues? How quickly time passes 
when we are having fun!” That is surely the secret of teaching mathematics - make it fun. 

I have taken the occasion to review the previous issues, and I am awed by the total 
volume of great unique and still-relevant ideas in those issues. They are collector’s treasures, 
and I cherish each and every one. 


We all look forward to each issue with great anticipation. It is astounding to learn 
such ingenious and effective ways that others have used this technology. 


So this toast of appreciation is not only to Josef and Noor, but also to all of the past 
and future authors. Your ideas are enriching mathematics for students, educators, and 
mathematical hobbyists all over the world. 


Aloha and Mahalo, 

David Stoutemyer 
Chairman of the Board 
Soft Warehouse, Inc. 



The 25th issue of the DERIVE Newsletter. 

This is a wonderful moment: The 25th DNL. This means that I have reveived 24 issues before. There 
are 4 issues a year. 24 divided by 4 gives ...wait a moment ... I can do this without DERIVE ... just 
another second ... 6 ... I think this gives 6. What were we talking about? Oh yes. 6 years. Yes, it is 
tme: The DERIVE User Group serves a growing community of DERIVE users and DERIVE enthusi- 
asts for more than 6 years now. 

The DNL is the major vehicle for communication among DUG members. And it serves its purpose 
very well. It is informative, challenging, and critical. The DNL is a valuable source of hints for the 
newcomers, it is a highly appreciated source of know-how for the experts, and it is a very serious and 
important means of feedback for the authors of DERIVE. I have seen both David Stoutemyer and Al- 
bert Rich sitting in their office in Honolulu studying the newest issue of the DNL very carefully - 
sometimes with deep wrinkles on their foreheads. 

I cannot remember a single one of the 24 issues having been published late. Knowing Josef personally 
for more than 8 years now this is no surprise. He is one of the most diligent and reliable persons I ever 
met. 

Josef and his wonderful wife Noor are absolutely charming people. If you haven’t met them yet per- 
sonally, come to one of the various international conferences (Josef lists them in the DNL) and you 
will make friends with them very quickly. There is no way not to become very friendly with them. 

DERIVE and the DERIVE User Group. The first would not be the same without the second. 

Josef and Noor, thank you! 

Bernhard Kutzler 

Managing Director of Soft Warehouse Europe 
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"A Day in the Life " 

DERIVE as a Demonstration Tool 
in Upper Secondary Mathematics 

Neil Bibby, Lancaster, UK 

During recent years I have developed an integrated use of the DERIVE package into 
my work at all levels of upper secondary mathematics. This has been primarily as a 
demonstration tool, by which I mean that the students have not had direct hands- 
on use of the package. There are teo good reasons for this: 

(1) DERIVE is a sophisticated piece of software which requires a thorough knowl- 
edge of mathematics on the part of the user, including a strong grasp of notation 
and notational variants. Although there is some evidence that it has been used suc- 
cessfully as a hands-on tool for students, equally other work has cast doubt on this. 
I think of DERIVE as a professional tool which in professional hands can provide a 
rich learning environment. 

(2) The availability of graphic calulators enables students to have their own graph- 
ics system in the palm of their hand: so my preferred way of working was to inte- 
grate the use of the two where this was possible. The students can create their own 
version of the DERIVE screen display to reinforce what they have seen on the 
DERIVE display. 

It will become apparent that my examples for the most part use the graphic capa- 
bilities of DERIVE. I do not claim that DERIVE is uniquely suited to this: no doubt 
much of what I present is equally implementable on other graphics systems. How- 
ever the graphics system of DERIVE is flexible and easy to use once one has mas- 
tered a few basic techniques. 


Example 1: n and Van Ceulen's algorithm 

(British Year 9 = Grade 8) 

It is surprising that the treatment of the topic of n is so little changed by the avail- 
ability of microtechnology. Apart from encounters with recurring decimals, it is the 
only place in elementary mathematics where the infinite process has to be really 
faced up to. Now that we have the computational tools there seems little to stop us 
from doing just that. The Archimedian approach of starting with a regular hexagon 
of unit radius and then doubling the number of sides is now eminently feasible on 
DERIVE. It's easy to draw regular n-gons as the following diagram shows. 



y 

1.5 

0.5 • V\ 

M x 

-2 -1 .S 

0.5 A 1.5 Z 

■ -0.5 ■ / / ■ 

-1.5 

COMMAND: Center Delete Help Moue Options Plot Quit Range Scale Transfer 

Uindou aXes Zoom 

Enter option Deriue XM 

Cross x:l y:l Scale x:0.5 y:0.5 ZD-plot 
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The screen shots above are from DERIVE for DOS from 1997. 


In fact the left diagram shows a hexagon and a 360-gonH The zoom facility of 
DERIVE allows us to easily see this in the right diagram. 

But how can we calculate n? Ludolphus van Ceulen (born late 16th century in The 
Netherlands) devised a simple iterative method based on the Archimedian approach 
- to this day n continues to be known as the 'Ludolphine number - die Ludolph'sche 
Zahl’ in German speaking Europe! By applying Pythagoras to two triangles in the 
following figure, we can simply relate the length s j + , of the side of a regular 2n-gon 

to the length s \ of the side of a regular n- gon. The relation is that 

s t +1 = ^2 - -^4 - s ; 2 . Starting with a hexagon of unit side, i.e. n = 6 and So = 1, we 

can approximate n as the semi-perimeter of the polygon by 3-2' ■ ,s t . The ITERATES 
and ITERATE functions of DERIVE make this very easy to do, and we can easily ad- 
just the accuracy of the calculation to avoid cumulative rounding errors occuring. 



f 2 k - 1 j 

#11: VECTOR (ITERATES(7(2 - 7(4 - s )), s, 1, 10)) -3-2 , k, 1, 11 

l k 

#12: [3, 3.1, 3.13, 3.13, 3.14, 3.14, 3.14, 3.14, 3.14, 3.14, 3.14] 

#13: PrecisionDigits 20 

#14: NotationDigits := 20 

f 2 k - 1 j 

#15: VECTORl (ITERATES(7(2 - 7(4 - s )), s, 1, 15)) -3-2 , k, 1, 16l 

#16: [3, 3.1058285412302491481, 3.1326286132812381971, 3.1393502030468672071, 

3.1410319508905096377, 3.1414524722854620752, 3.1415576079118576454, 3.141583892148318406, 
3.141590463228050057, 3.1415921059992713602, 3.1415925166921566904, 3.1415926193653861406, 
3.1415926450336881605, 3.1415926514508704393, 3.1415926530552116196, 3.1415926534558914124] 


#17: n 


#18: 3.1415926535897932384 
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Example 2: Population growth and exponential modelling 

(British Year 11 = Grade 10) 

This work formed the basis of some GCSE coursework. We started with the popula- 
tion data for the major Italian cities as given in the following table: 


Coioa 

Turin 

Milan 

Rome 

Naples 

Palermo 


Tablet Population of the main Italian cits 


rSoo/i 

i8jo/t 

t88o/i 

1870/1 

1880/I 

1900/1 

1910/t 

V 

110 

x *9 

130 

180 

*35 

* 7 * 

7 « 

135 

17S 

108 

*54 

33 * 

4*7 

*35 

* 4 * 

* 4 * 

161 

3 ** 

4*3 

579 

1*3 

*75 

1S4 

*44 

300 

4*3 

54 * 

4*7 

44 ? 

4x7 

449 

494 

5*4 

7*3 

1}> 

I So 

1 u 

1X9 

*45 

3X0 

34 * 


s (m thousands) itoo-lfSo 


910/1 

I 9 JO/I 

1940/ 1 

1950/ 1 

1980/1 

1970/1 

1980/1 


808 

*35 

*4* 

784 

8ti 

7 *o 

<01 

597 

**9 

711 

10*8 

ri78 

1104 

* 3 * 

99* 

1 1 16 

1*60 

ij«j 

17*4 

t«jj 

* 9 * 

1008 

1156 

185 1 

tilS 

1800 

1831 

/** 

*39 

166 

totx 

1X83 

**33 

HEX 

3*4 

390 

4x1 

49 X 

588 

851 

700 


Noce the continued growth of dries in the inter-war years (despite fascist attempts to encourage ‘ruralityl and the surge in the 
period 1950-70. Naples was Italy's largest dry until the First World War. 

Source: B. R. Mitchell. International Historical Statistics, Europe njo-tpit (N.Y.. 1991). 


We decided to concentrate on Milan from 1880 to 1980 (not necessarily the easiest 
choice - why?) the data can be plotted as follows: 



1800 

population 

1600 

1400 

1200 

1000 

800 

600 m 
400 
200 

(time - 1800)/ 10 


123456789 10 


The problem now was to try to fit a curve of the form y = a- b x to this data, with y 
as the population in thousands and x as the number of decades from 1880. Initial 
suggestions usually were a = 322 and for the choice of b to give the visual best fit to 
the other points. Since the doubling- time appeared from the table to be about 4 
decades, by the rule-of-70 this gave a growth-rate of about 18% (since 70/4 « 18), 
and hence b « 1.18. 

Further discussion led to the method of visual-best fit being questioned, and the TI- 
82's statistical ExpReg function being used to provide a least squares fit. I hasten 
to add that the intricacies of the least-squares method were not discussed: the TI- 
82 function was therefore treated as reliable and trusted ’blackbox’ which would 
give an appropriate result. (DERIVE's FIT function could also have been used for 
this purpose). The values obtained were a = 412.5, b = 1.198 (4 s.f.) and plotting 
y =a-b x with these values gave the following result: 
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0 

LN(322) 

> 


1 

LN(493) 



2 

LN(579) 



3 

LN(836) 



4 

LN(992) 


[x, a-x + b] , 

5 

LN(1116) 



6 

LN(1260) 



7 

LN(1583) 



8 

LN(1724) 



9 

LN(1635) 

> 


#26: 0. 180389 -X + 6.02224 

0 . 180389 *X + 6.02224 

#27: e 


The visual immediacy of this treatment, combined with the use of TI-82 (albeit in a 
'blackbox' mode) proved pedagogically. 


Example 3: Simpson's Rule 

(British Year 12 = Grade 1 1) 

My preference with this algorithm is to emphasise the fitting of a quadratic function 
by three points of the graph of the function whose integral is required. A nice exam- 
ple to start with is INT(cos(x), x, -n/3, ji/3). The interpolating quadratic is easily 

9x 2 

seen to be v = 1 . 

2 k 1 

The graph of this gives the following, together with that of cosine; omitting here but 
also worth including, is the graph of the second degree Taylor approximation of co- 
sine at 0, to emphasise that this approximates better at 0, but doesn't interpolate at 
all. 



A comparison of the corre- 
sponding integrals is then illu- 
minating: the exact value is V3 
(1.732...), the Simpson value is 
1.745.... and the ’Taylor’ value 
is 1.711.... 


Students should appreciate that Simpson's Rule (actually due to Newton originally - 
who else?) is exact for quadratics: this seems obvious, but a few corroborative cal- 
culations using DERIVE push home the point. However one of the most intriguing 
features of Simpson's Rule is that it is exact for cubics!! This 'two-for-the-price-of- 
one' aspect needs some explanation. A formal proof is possible of course, but an 
informal graphical approach might be appropriate and convincing initially. 

The diagram shows the graphs of a cubic, its quadratic interpolator, and the differ- 
ence of the two. This difference is a new cubic, whose zeros must be a, a + h, and b 
(= a + 2h). 
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So the point where x= a+ h must be a point of inflexion, i.e. the symmetry point of 
the graph. Hence the two 'lobes' of this cubic are of equal area, and thus Simpson's 
Rule is exact!! The visual plausibility of this is compelling, but should not detract 
from the details of the argument; once again the two aspects corroborate each 
other. 



Example 4: Polar and Cartesian Coordinates 

(British Year 13 = Grade 12) 

DERIVE version 3 has the facility to plot implicitly. This is naturally slower than the 
explicit plotting, but can be used to good effect, especially when explicit plotting is 
simply not possible (e.g. y 5 + y = x). It is especially effective when the equivalence 
of polar forms and cartesian (rectangular) forms of an equation are being consid- 
ered. 


One approach to the conic sections is to use the focus-directrix definitions to arrive 

atr = . (It is instructive to consider the limagons r = a{\ + e-cost? ) 

1 + e-cos£? 

first. As an example we consider r = , where a = e = 0.5. This gives an el- 

2 + cost? 

lipse (or hypObola, as I sometimes call them!), and is easily plotted using DERIVE: 
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It is now instructive to ask "what is the cartesian equation of this curve?" We use 
the substitutions r = + y I 2 and cos 0 = ——^= to obtain the cartesian equiva- 

v* 2 + y 2 

lent f(x + }) 2 +3 y 2 =1. The static graph above cannot do justice to the dynamic effect 

of seeing this plotted (remember to switch back to rectangulars!) in a new colour 
exactly on top of the polar verson. The monochrome calculator displays cannot 
compete here, even if they could plot implicitly! The cartesian form provides us with 
new information: we can see from the equation that x = -j is a line of symmetry, 

and that the lines y-±-j^ are tangent lines: in the follwing diagram we confirm 
these: 


2 2 1 

7(x + y ) = 

x 

# 60 : 2 + 

2 2 
7(x + y ) 

2 2 

2 2 7(x + y ) 

#61: 7(x + y ) = 

2 2 




2-V(x + y ) + x 







r 

2 2 ^ 







2 2 

</(x + y ) 

2 2 





#62: 

II 

>> 

+ 

X 


.(2.V(x + y ) + x) 


2 

2 2 

2 2 2 



2 2 


#67: 

4- (x 

+ y ) 

>, 

+ 

X 

/-\ 

i—i 

I 

X 

II 


w 

2-V(x + y ) + x J 











2 

2 2 

2 2 2 


2 2 

2 2 2 

! 2 


4- (x 

+ y ) 

= (x - 1) -(x + y ) 

#63: 

X-V(X + y ) 

+ 2-x + 2-y = V(x 

+ y ) 

#68: 











2 2 


2 2 

2 2 

2 2 2 2 




x + y 

#64: i 

(x -7(x + y ) 

+ 2-x + 2-y = 

: + y )) - *./(x + y ) 










2 

2 

2 


2 2 

2 2 


#69: 

4- (x 

+ y ) = 

(x - 1) 

#65: ; 

2-x + 2-y = 

(1 - x)-/(x + y ) 











9 f 

1 1 

2 2 


2 2 

2 2 

2 

#70: 

— * 

+ - 

+ 3-y =1 

#66: i 

(2-x + 2-y 

= (1 - xW(x + y )) 


4 l 

3 J 




I hope that these examples will convey the favour of this approach. I now find 
DERIVE an indispensible tool in my day-to-day teaching. 

Please don't hesitate to contact me if you need further advice or assistance. 

Neil Bibby 

Department of Mathematics 
University College of St Martin 
Lancaster LAI 3JD, UK 
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DCRITflflf 

TANGRAMS with DERIVE 

Alfonso J. Poblacion, Valladolid, Spain 

INTRODUCTION 

One of the oldest parts of Recreational Mathematics concerns Dissection problems. Among them, 
the ones derived from seven pieces (called tans) known as tangrams are very popular, (see figure 1). 
Arranging these tans it is possible to form a wide variety of figures. The rules are simple: the seven 
tans must be used and it is not allowed to overlap any of them, although we can let holes among the 
pieces. Martin Gardener in [1] classifies in three types the questions we can deal with tangrams: 

1 . to find how to build a given tangram or to prove its impossibility 

2. to find how to form several real figures such as animals, objects or persons in the most artis- 
tic or amazing way 

3. to solve problems about Combinatorial Geometry that the seven tans set forth. 

In the mentioned reference we can find more about these categories of problems and the history 
of tangrams. 


HOW TO USE DERIVE 


We will try to play with the tans using DERIVE and our imagination. First of all, we will define 
the seven pieces (lines#3 and #4) from its original place in the initial square. 


# 1 : 

# 2 : 


#3: 


#4: 


[Notation := Decimal, NotationDigi ts := 3] 


The seven tans 




0 

0 


0 

2-72 ' 


J2 

72 


72 

72 

a := 

0 

2-72 

, b := 

2-72 

2-72 


0 

0 


0 

2-72 






3.72 

72 





2 

2 


' 72 

0 


2-72 

72 

d := 

2-72 

0 

P ’ “ 

3-72 

3.72 


2-72 

J2 

P •" 

2 

2 


. 72 

0 


72 

72 





3.72 

72 





2 

2 


0 

0 


72 

72 

72 

0 


2 

2 

3-72 

J2 


3-72 

72 

2 

2 

, g := 

2 

2 

72 

72 


72 

72 

2 

2 


72 

72 

0 

0 


2 

2 



D-N-L#25 


AC DC Tour 


p 45 



You might wonder about the colours with the 
DERIVE plot. This is easy to perform with 
DERIVE 6. 

Transfer by applying Mark and Copy the 
figure into the Algebra Window. Mark the 
graph and activate “Convert Picture Object” 
in the Edit-Menu. Then you can switch to 
Paintbrush and treat the graph as you like. 

(See also Tania Koller's note in DNL#63, 

0.5 l 1.5 2 2.5 3 page 14). 


#5: Two basic functions to play 

#6: T(x, y, «, u, v) := [COS(o!)-x - SIN(«)-y + u, SIN(o<)-x + COS(«)-y + v] 

MOVEMENTfw, «, u, v) := VECTOR(T(w , w , a, u, v), i, DIMENSION(w)) 

#7 : i , 1 i , 2 

The first one moves a point (x.y) by a rotation of angle a followed by a translation of vector 
(u,v). The last applies a movement T to a whole tan w. As I said before we can add REFLECTIONS 
or any other affine motions we desire. 

Now the goal is to find the movements needed to bring the tans from its initial coordinates to the 
new ones required to form a figure. Let me explain it with an example. We desire to sketch the Yacht 
in figure 2. After discovering where the tans must be put (this can also be tried using DERIVE ), but it 
can make us spend a lot of time), we can find the movements by testing or calculate them exactly. I 
got the Yacht with these. (Simplify and plot the result in Connected Mode): 



p 46 


AC DC Tour 


D-N-L#25 


#13: MOVEMENT 


3 -n 


b, 


12 


2-V2-COS 


f 1 1 

f 1 1 

5 

r~ 

i— i 

— 

+ 2-^2 -SIN 

+ , 2* ,/2. COS 


1 12 J 

l 12 J 

2 

, 12 > 


2-J2 -SIN 


(”)) 


#14: MOVEMENT 


d, ATAN 

2 


H“hr) * S “(V)) 


r r \ \ 


SIN 


^ 6 


* IHt) - D 1 




4 

’ 4 SIB (t) - 


k 6 


+ 1 . COS 


COS 


L 3 


4-(« 

■) 


f 12 


- SIN 


(V)] 


j ( sim (t) * ‘) 


5 r o 

+ , (2 -J2 - 2). COS 

2 l 12 


f 1 ^ 

• SIN 

l 12 J 



Of course the last ones cannot be obtained at first sight. (The complete procedure can be found in 
file DERITANS.MTH.) 

The next two functions can be helpful to find angles and distances: 


x_ *(y_ - x) - y_ • (x_ - x) 
2 1 2 1 


LINE(x_, y_) := 


y_ - x_ 

1 1 


DIST(x_, y_) := V 


(y_ - x_ ) + (y_ - x. 

11 2 


v) 
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You can find much more about tans in the mentioned reference. To practise I propose you two 
different problems: 

(1) Try to find the movements to compose these figures: 



(2) What is the biggest area that a tangram can contain? (In [1], it is given a solution but as far as I 
know it is not still proved that this is the biggest one. This is an open problem). 



Finally, two crazy suggestions. It could be funny to play with 
DERIVE and the tans among several players trying to find the 
movements before the opponents. Or perhaps anyone could think 
about setting out exams in which the students are given a figure 
to find the movements needed to be formed. So our students 
would love and enjoy mathematics much more or perhaps they 
hate us even more. If someone dares to put it in practice, tell us 
how it goes. 


REFERENCES 


[1] Gardner, M. Viajes por el tiempo otras perplejidades matematicas. Capitulos 3 y 4. Ed. Labor. 
S.A., 1988 (English version: Time travel and other mathematical bewilderments). 


I was really fascinated by Alfonso 's article. It 
combines fun and creative phantasy with 
mathematics in a wonderful way. I imagined 
wooden tans to play with. So I tried to pro- 
duce shaded tans - just for fun - to present 
the solutions or the problems in a nice form. 
That is the "wooden" Yacht. If you want to 
work with the shaded tans then use the 
new tans. 

(The complete solution for the Yacht will be 
found on the diskette of the year. I hope to 
be able to present the actual files to be 
downloaded on the homepage of my school 
in the near future.) 

Josef 
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A1 Rich's interesting answer to Nigel Backhouse's comparison of calculation times: 

The following is in response to Nigel Backhouse comments dated 27 January 1997 to the DERIVE News mail- 
ing list concerning integer factoring of 2 A 101-1: 

To factor integers, DERIVE uses trial division, followed by a Monte-Carlo Primality test, followed by the Pol- 
lard Rho algorithm as detailed in Knuth, The Art of Computer Programming, Volume 2. For the most part, 
these algorithms are implemented using loops rather than recursion. This combination of factoring methods is 
usually the most efficient until the size of the second largest prime factor exceeds about 12 digits. 

Glossing over details, the more complicated elliptic curve method is then usually the most efficient until the size 
of the second largest prime factor exceeds about 25 digits. Quadratic residue methods are usually the most effi- 
cient beyond that. 

The times Backhouse reports measure the relative speed of the algorithms used by the various computer algebra 
systems to factor 2 A 101-1, more than they measure the speed of the underlying arithmetic. To try to determine 
the underlying combination of algorithms being used by the various systems, construct a family of integers hav- 
ing n digits in its second largest prime, and then compare the times required to factor the family of integers by 
the various systems. 

The relative rankings on different CPUs are likely to be the same, and more or less proportional to sequential 
memory fetch and store, so testing on additional CPUs is not likely to yield much additional information. 

The bignum arithmetic in DERIVE uses the same nonrecursive algorithms (also described by Knuth) as your 
other systems, but DERIVE's is implemented in Intel 32-bit assembly language. Therefore, the DERIVE big- 
num arithmetic should be faster than bignum arithmetic compiled from languages such as C. To test this, com- 
pare times for multiplying some non-special numbers having 100, 200, 400, ... digits, and time the multiplica- 
tion. Times for the smallest examples are likely to masked by the timer resolution unless done using ITERATE. 



Here they are - at the occasion of Josefs and Josefs daughter’s birthday (same day! !): 
1 st row: Constantin, Dominic, Yvonne, Kim (the baby from page 1) 

2 nd row: Moritz, Naomi and Maxime 
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From the picture we can see that your friend will be in position to catch the ball when t = n/2. We will 
draw your friend's path from the time the ferris wheel starts until the ball is supposed to arrive. 

In DERIVE author: [20 sin(t/ 2 ), 25 - 20 cos(t/2)] 

To plot this enter the Plot Window and choose the Plot option. DERIVE will prompt you for the range 
of the parameter, t in this case, enter 0 for Min and jc/ 2 for Max. An eighth of a circle will be drawn. 

What about the ball? Let P(t) and Y(t) be the horizontal and vertical components of the position 
of the ball. Let 6 be the angle that the ball is pitched (angle with the negative x-axis). Then 

P(t) = 28-64 cos(<9) t 

Y(t) = 5 + 64 sin(#) t - 16 1 2 

(A nice question for the non-English or non-American students: Where does the number 1 6 come 
from?? Josef) 

We need to determine the choice of <9 that will cause the curve determined by [P(t), Y(t )] to intersect 
the path of our friend when t = n/2. We will determine the time delay later. Author 

[28 - 64-COS(0)-t, 5 + 64-SIN(0)-t - 16-t A 2] 

In DERIVE (for DOS) 9ist typed as Alt-h. We can not draw any graphs at this point. We must make a 
choice for either 9 or t. This is the beauty of having a tool like DERIVEl We can experiment. Let's 
make a naive choice, 9 = tt/ 4. Use Manage/Substitute to draw, replace 9 and then Plot the resulting 
expression for t going from 0 to 1 . As you see on the graph below, this overshoots the mark and is too 
long a time. Now go back to the Algebra Window and use Manage/Substitute to replace 9 in our 
original expression with 71 / 8 . This time when you Plot let t run from 0 to 0.7. As we see this under- 
shoots the mark and is still too long a time. 



You continue the investigation and find out the value for 9 that will work for this problem. Hint: 7i/8 
isn't off by much. Once you have the correct angle, you can use the trace to determine how long it 
takes the ball to reach your friend's path. Now you can decide how long to wait. If you want DERIVE 
to solve for t analytically, you can also do that once you have the value for 6. Go back to the Algebra 
Window and substitute n/2 into the parametric equations for your friend's position and then set P(t) 
with 6 given as the value you found from your graphical exploration equal to the x-coordinate of your 
friend's position. 

This example shows that it is possible for students with only a knowledge of trigonometric rela- 
tions to solve a very sophisticated problem using that knowledge and some good graphics. We grant 
that some analytical steps were bypassed and the answer may be an approximation to the ’’correct” 
answer, but was the student any less involved with the mathematics? We leave that question for your 
decide. 
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The TI-92 Corner 
(Edited by B. Waits, F. Demana, B. Kut zl er & J. Bohm) 


DERIVE and TI-92 do not always yield the complete set of solutions to 
systems of linear equations containing parameters 

Karl-Heinz Keunecke, Kiel, Germany (kh@KeuKiel.NetzService.de) 

At school in the field of linear algebra it is often helpful to make use of CAS to solve systems of linear 
equations. 

After the input of the augmented matrix students can calculate the row echelon form of the matrix by 
means of the DERIVE command Row_Reduce or with the command rref() of the TI-92 from which 
the solution of the linear equation can be read directly. This is shown in the following example. 

2x + y-3z = 4 i a a\ nncA x = 2 


2 1 

-3 4 ) 

l 

0 

0 

2 ) 

1 -1 

1 0 => 

0 

1 

0 

3 

3 -1 

-1 2 

0 

0 

1 

1 


x + y + z = 0 => 1-1 1 0 => 0 1 0 3 v = 3 

a 0 l 3 — 1 -1 2 lo 0 1 lJ 1 

3x-y-z = 2 k J k J z = 1 

This method can be applied to equations with parameters and is shown in the example below in which 
t is an arbitrary real number. 

3x-2y + z-2t 3-2 1 2t T fl 0 0 t + l) x = t + l 

5x-4y-z = 2 =^> 5 -4 -1 2 => 0 1 0 t + 1 =^> y = t + \ 

z — t- 1 


x - 3y - 2z = 2/ + 6 


(3 

-2 

i 

2 1 

1 1 

0 

0 

t + l) 


5 

-4 - 

-l 

2 

=> 0 

1 

0 

t + l 


u 

3 - 

-2 

21 + 6 

J lo 

0 

1 

t-l) 


F£ 

■■ghIO 

T“ 

Iciear 

i 

a-; 

n 

However, 

if 

one 

selects 

the 


2 b "3 
>43 3 

.4 -3 9 

IT 2 6 

' rref 4 3 


LG G G 
'2 6 "3 

4 3 3 

.4 -3 9 

' 1 G 3/; 
0 1 "I 

.GGG 


2x + 6y - 3z = -6 
4x + 3y + 3z = 6 . 
4x-3y + 9z = £ 


'1 G 3^2 3 ] 
Gl-l -2 

.0 0 0 0 . 


U.4 "3 

9 k JJ 

T2 6 

-3 -6T 

4 3 

3 6 

J.4 "3 

9 k J, 


iTTWi? — — y~f 5 ^y fh-t y fe ~y Ffi then neither DERIVE nor the TI-92 give the correct 

— k j] [g Q Q ^ (complete) solution. In the display the command rref() 

[T 2 6 " 3 ~ 6 T| was applied to the augmented matrix of the system of 

■ rref 4 3 3 6 | k = IS 

U.4 -3 9 k JJ equations. From the result it can be concluded that for 

g l ^ l 2 ^-2 ^ E R e the system of equations has no solution. The 

■-■■■■■ _, Q , bottom row may be interpreted as 0x+0y+0z=l and 

irniN rhomis ruHciB/30 gives accordingly a false statement. 

This result is obviously incorrect. Because for when £=18 you get a different solution as can be seen 
from the adjacent display of the screen. The solution is read as, 

3 

£ = 18: x = - — z + 3 A T = z- 2, ze Re, where z is an arbitrary parameter. 

The question is then, why the solution for £ = 18 is not recognized by the CAS. The cause is, of 
course, the algorithm of the respective CAS. But even without knowing this, one can understand why 
the above solution is not recognized. For this one merely needs to determine the solution step by step 
by means of elementary row manipulations. This can be carried out quite easily with the commands of 
the TI-92. By means of the commands mRow (factor, matrix , n ) one can multiply the row n of matrix 
by factor. The command mRow Add {factor, matrix, nj,n 2 ) allows the row nj of matrix to be multiplied 
by factor and to be added to the row n 2 . 
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1.4 "3 

1 nRowRddl 


lear a- z... 


9 

n 

L4 -3 


[2 6 -3 -61 

■ 

"2, 

4 3 3 6 

,1,2 


.4-3 9 k . 



6 

-9 

-3 


-6' 
IS 
k . 


,-G][4,3,3,G][4,-3,9,k]],l,2> 


FUNC E/30 


1.4 -3 9 k J| 


Ffi 

ear a-z... 


1 nRowfldd 


"2, 


6 

-9 

-3 


-3 -6 

9 IS 


=■ 1,3 


k 

2 6 -3-6 

0 "9 9 IS 

G -15 15 k+12. 


..][□, "9,9, 18][4, “3,9, k]],l,3> 


RAD AUTD 


FUNC 3/30 


Firstly row 1 is multiplied by -2 and is then 
added to row #2. 

Then row #1 is multiplied again by -2 and the 
added to row #3 (below left). 

Now row #2 is di-vided by -9 (below right). 


ifl^T FT^ Y F]:t Y FHt Y F? Y Ffi 
K fr— |R1 gebra|Calc|Qther|FrgnIQ|Cl ear a-z... 


LG "15 15 k+12, 




'2 

6 

-3 

-6 


■ nRow 

- 1/9, 

G 

-9 

9 

IS 

,2 



.G 

-15 

15 

k + 12. 




'2 6 -3 -6 

G 1 -1-2 

.0 -15 15 k+12. 



L 


MAIN RAD AUTD FUNC H/30 



iTiT^ril Y F:=t Y FHt Y FE Y Ffi 1 

f— |Rlgebra|Calc|Ot.her|PrgnIO|Clear a-z... | 




Lg 

G 

G 

k - 1SJ 



r 1 G 3/2 

3 




■ nRow 

1 

k - IS ? 

G 1 -1 

-2 


? 3 




.0 G G 

k - 

is. 







ri 

G 

3/2 3 ] 





G 

1 

-1 -2 





.0 

G 

0 1 . 

HAlDjUUliMAiDjlDjlBA 

MAIN 

RAD AUTD 

FUNC B/30 



If row #2 is multiplied by 15 and then added to 
row #3 the result shown on the left is obtained. 
The last element of the bottom row is £ - 18. 
Therefore it is necessary to differentiate between 
the two cases £=18 or £^18. 

In the first case the set (1) of equations is solu- 
ble but not in the second. The solution for the 
case £ = 18 can, after further transformations, be 
read directly from the display: 

3 

£ = 18: x = — — z + 3 a y = z-2, zeRe 
2 

The CAS have not made the distinction between 
the cases, but rather have erroneously divided the 
last row by £ - 18. Even with the TI commands 
this is consistently possible as can be seen in the 
display on the left. This result is identical to the 
result of the command rref() at the beginning of 
this report. 


The lack of distinction in cases and illegal division respectively can occur only, if in the row echelon 
form of the augmented coefficient matrix, in any single row only zeroes and, to the far right, an ex- 
pression appear. The set of solutions is then either empty or it has a dimension greater than or equal to 
1, as if shown in example (1). 

Remarks: 

1. The reported problem arose in the course of vector algebra. The students should find out how 
the intersection of the planes 

Ei: 2x + 6y-3z = -6 

E2: 4x + 3y + 3z = 6 depends on the real number £. 

E 3 : 4x-3y + 9z = £ 

2. The commands mRow() und mRowAdd() are very well suited to helping students master the 
method of solving systems of equations. They can concentrate fully upon the method, since the 
calculator relieves them of the numerical work. 
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The T I - 9 2 Corner: Script & Program 


D-N-L#25 


Example of a script (in combination with W. Propper's program ABLEIT() - find it among the 
files!) 

: We investigate the differential quotient of a special function for x=0 
:We set up the screen: 

C : s e t mo d e ({"Split Sc r een" , " Top- Bot t om" , " Spl i t 1 App" , " Home" , " Spl i t 2 
A p p " , "Text E d i tor", "Spl i t Screen Rati o " , " 2 : 1 " } ) : cl r h o me 
C: ( abs ( x) / 2- 2) A 2»f ( x) 

: f ( x =0 ) = ? 

C: f ( 0) 

:we define the rate of change approaching from right hand side: fromr(x) 

C: ( f ( x+h) - f ( x) ) / h»f r omr ( x) 

: Now special for x =0 : 

C: f r omr ( 0) 

: F r o m r i g h t : I i me s mi t h » 0 
C: I i mi t ( f r omr ( 0) , h, 0) 

: Wh y that? 

:h mus be positive! Hence: 

C: f r omr ( 0) | h>0 

: Fai I ed ! Hel p the Tl and try Expand 
C:expand(fromr(0))|h>0 
: t h e I i mi t is? 

C: I i mi t ( a ns ( 1) , h, 0) 

: t he slope of the tangent is? 

:the intersection of the tangent a 
: t he equat i on of t he Tangent ? 

C: 5 2 x + 4 » 1 1 
: I et us plot it! 

C: cl rgraph 

:the settings for the Graph screen 
C : 5 1 2 » x mi n : 1 2 » x ma x : 3 » y mi n : 7 » y ma x 
C: gr aph 1 1 

:Now let's approach from the left: 

C : ( f ( x ) ■ f ( x - h ) ) / h » f r o ml ( x ) 

C:expand(fromi(0))|h>0 
C: I i mi t ( a ns ( 1) , h, 0) 

:the derivatives are different -> not differentiable!! 

: we obtain another tangent 
C: 2 x + 4 » 1 2: graph 1 2 

: C o u I d you imagine the curve's shape? 

C: gr aph f ( x) 

: I t has a "Cusp" 

: Do you have an idea how the derivative looks like? 

: Zeros of the derivative? 

: What's the value of the derivative for x=0? 

: "Let" calculate the derivative and store under dqf(x) 

C : 1i ( f ( x ) , x ) » d q f ( x ) 

C: gr aph dqf ( x) 

C: dqf ( x) | x<0 
C: dqf(x) j x>0 

: Can you explain the last results? 

: Do n t forget house keeping!! 

C : del v a r dqf,fromr,froml,dqf,tl,t2 
C: set mode ( "Spl i t Screen", "Ful I ") 

: Back wi t h Di amond HOME! 

I used this teaching unit to repeat the differentiability of a function using a View Screen. The two prod- 
ucts - the script and the program - supported each another in an excellent way. Josef 










